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ABSTRACT 

New observations show that the hghtcurve of Kuiper belt contact binary 
(139775) 2001 QG298 has changed substantially since the first observations in 



m 

^' 

p I 2003. The 2010 lightcurve has a peak-to-peak photometric of range Am2oio = 

Q 0.7 ± 0.1 mag, significantly lower than in 2003, Am2oo3 = 1-14 ± 0.04 mag. 

-(— > This change is most simply interpreted if 2001 QG298 has an obliquity near 90°. 

c3 The observed decrease in Am is caused by a change in viewing geometry, from 

equator-on in 2003 to nearly 16° (the orbital angular distance covered by the 

^ object between the observations) off the equator in 2010. The 2003 and 2010 

C lightcurves have the same rotation period and appear in phase when shifted by 

1/-J an integer number of full rotations, also consistent with high obliquity. Based 

^ on the new 2010 lightcurve data, we find that 2001 QG298 has an obliquity e = 

L^ 90° ± 30°. Current estimates of the intrinsic fraction of contact binaries in the 

^-H Kuiper belt are debiased assuming that these objects have randomly oriented 

spins. If, as 2001 QG298, KBO contact binaries tend to have large obliquities, a 

larger correction is required. As a result, the abundance of contact binaries may 

be larger than previously believed. 

Subject headings: Kuiper belt: general — Kuiper belt objects: individual (2001 
QG298) — methods: data analysis — techniques: photometric 



1. Introduction 

Located approximately between 30 and 50 AU from the Sun, the Kuiper belt represents 
the outer frontier of the currently observable solar system. The belt is estimated to hold 



roughly forty thousand objects larger than 100 km in diameter ( Jewitt fc Luu|[l995 Trujillo 
et al.||2001 1 [Bernstein et al.|2004[[Fuentes fc Holman|2008[ [Eraser et al.||2008D . More than one 
thousand objects have been detected but only six hundred have multi-opposition observations 
and well-determined orbits. The Kuiper belt objects (KBOs) are remnants of outer solar 
system planetesimals and their study may shed light on the complex process that led to 
the formation of planets. Our understanding of the physical properties of KBOs remains 
rudimentary as most are too faint for detailed investigation. 

Many KBOs are binary. Among the dynamically cold, "classical" population (low ec- 
centricity and inclination KBOs located between the 3:2 and the 2:1 mean-motion resonances 
with Neptune, at 39.4 and 47.8 AU) the binary frequency is inferred to be ~20% whereas in 



other dynamical subsets the fraction is lower, at 5 to 10% (Stephens &: Noll||2006 Noll et al. 



2008). Binaries probably formed early on (Weidenschilling 2002), before the Kuiper belt 



lost more than 99% of its original mass (Kenyon & Luu 1998 Booth et al. 2009), because 



the current-day number density of KBOs is too low for frequent pair encounters. It is also 
possible that planetesimals formed in clusters of two or more bodies and that the binaries 



we see today are direct remnants of that process (Nesvorny et al. 2010). Thus, binaries are 



particularly useful probes of the conditions under which planetesimals formed. For instance, 
the symmetry in surface properties of KBO pairs argues for a chemically homogeneous for- 



mation environment (Benecchi et al. 2009). But the current binary abundance also bears 



record of the effects of coUisional and dynamical processes that have eroded the original 
population ( Petit fc Mousis|[2004 ). It is highly desirable to find ways to isolate the effects of 
formation, evolution, and erosion. The relative frequency of binaries within different regions 
of the Kuiper belt, and their distributions of orbital properties and relative mass of the com- 
ponents are quantities that should be sought to help clarify the formation and evolution of 



binaries (Schlichting & Sari 2008b Murray-Clay & Schlichting 2011 Nesvorny et al. 2011). 



Most known binaries are resolved (physical separations >1500 km; Noll et al. 2002 



Stephens &; Noll 2006; Stansberry et al. 2006 Grundy et al. 2011). The exception is the 



500 km X 175 km contact binary (139775) 2001 QG298, identified from analysis of its rotational 
lightcurve (Figure [l| [Sheppard fc Jewitt 2004, hereafter SJ04). Using time-resolved mea- 
surements taken mainly during 2003, the authors found that 2001 QG298 displayed extreme 
photometric variability (Am ~ 1.2 mag) and a relatively slow rotation period (P ~ 13.77 
hr). The large photometric range, caused by the eclipsing nature of the binary, exceeds 
the Am = 0.75 mag produced by two spheres in contact (axis ratio b/a = 0.5) because 
the contact binary components are tidally elongated along the line connecting the centers 
(Weidenschilling 1980 Leone et al. 1984). The large Am found for 2001 QG298 also implies 
that the system was observed almost equatorially in 2003. Models of 2001 QG298's lightcurve 



based on figures of hydrostatic equilibrium lend strong support to these assertions (Taka- 



hashi & Ip 2004 Lacerda & Jewitt 2007 Gnat & Sari 2010) and allow the bulk density of 



2OOIQG298 to be estimated at p = 0.59_q;q5 gem '*. This surprisingly low density implies 
that 2001 QG298 is mostly icy in composition and is significantly porous. 

The example of 2001 QG298 suggests that between 10% and 20% of KBOs could be 



contact binaries, depending on surface scattering properties (SJ04; Lacerda & Jewitt 2007). 
This high fraction is consistent with the observation that the frequency of resolved binaries 



increases dramatically with decreasing orbital separation (Kern & Elliot 2006). However 



despite their high abundance, contact binary KBOs remain elusive. The reason is that these 
objects are unresolved and can only be identified if they happen to be nearly equator-on. The 
consequence is that as many as 85% of contact binaries may go unnoticed due to unfavorable 
observing geometry. High phase angle observations that have been shown to increase the 



chance of detecting contact binaries (Lacerda 2008) can not be applied to the distant KBOs 



as their phase angles remain below 2° when observed from Earth. Interestingly, the Jupiter 
Trojans also contain a fairly large fraction of contact binaries, comparable to that found 



for KBOs. In addition to the well-known case of (624) Hektor (Hartmann & Cruikshank 



1978 Detal et al. 1994 Cruikshank et al. 2001 Lacerda & Jewitt 2007), two other contact 



binaries have been identified in a survey of 114 Trojans (Mann et al.||2007). When corrected 



by the probability of observing these objects close to equator-on, Mann et al. (2007) find an 



intrinsic contact-binary fraction of at least 6% to 10%. The apparently high abundance of 
contact binaries and the prospect of measuring their densities makes these primitive objects 
particularly interesting targets of study. 

QG298 has travelled a significant angular distance (~ 16°) in its heliocentric orbit since 
the first observations in 2003. In this paper we compare new 2010 lightcurve data with 
the SJ04 observations and make inferences on the obliquity of 2001 QG298- In Section 2 we 
describe our observations, in Section 3 we report our findings on the lightcurve and obliquity 
of 2001 QG298, in Section 4 we discuss the implications of our findings for the abundance of 
contact binaries, and in Section 5 we present our conclusions. 



2. Observations 

Visible-light observations were taken at the 2.5-m Isaac Newton Telescope (INT) oper- 
ating on the island of La Palma in the Canary Islands, Spain. The INT was equipped with 
the Wide Field Camera (WFC), mounted at the f/3.29 prime focus. The WFC consists of an 
array of 4 thinned EEV 4Kx2K charge-coupled devices (CCDs). The CCDs have an image 
scale of 0.33 arcsec per pixel and the full array covers approximately 34 arcmin squared on 
the sky. Our observations made use of CCD #4 which is the most central with respect to 



the field of view. Details of the observing times and conditions can be found in Table [T} 

Our target, 2001 QG298, was observed through the Sloan-Gunn r filter (Isaac Newton 
Group filter #214) which approximates the r band in the Sloan Digital Sky Survey (SDSS) 
photometric system. Consecutive 600 s exposures were obtained in a consistent manner to 
measure variations in the brightness of QG298- The rate of motion of 2001 QG298 across the 
sky was <2.32 arcsec/hr (0.39 arcsec in 10 minutes) keeping it well within the seeing disc 
during the observations. All frames were bias-subtracted and flat-fielded using the median 
of a set of dithered exposures taken during twilight. In the reduced images, the brightness 
variations of 2001 QG298 were measured differentially against 8 bright, not saturated, field 
stars to offer protection from fluctuations in atmospheric transmission and seeing. The 
mutual relative photometry of the comparison stars is stable to better than ±0.01 mag 
throughout all the observations. The field containing 2001 QG298 (R-A = 0^838, dec = 
1?964) at the time of our observations is included in the SDSS photometric catalog. This 
allowed us to absolutely calibrate our observations directly against SDSS magnitudes of 5 
of the 8 comparison stars (see Table [2]) chosen for their spatial distribution and exceptional 
photometric stability. 



3. Results 

3.1. The lightcurve of 2001 QG298 

We used our photometric measurements of 2001 QG298 to construct a lightcurve (mag- 
nitudes versus time) for this object. Figure [2] shows how our 2010 lightcurve and the SJ04 
lightcurve compare. Overall, the lightcurves have the same period and appear aligned in 
rotational phase. The 2010 photometric range is smaller by about 0.5 magnitudes. 

The apparent alignment in rotational phase between the 2003 and 2010 lightcurves is 
real. All times were measured relative to 2010 Aug 13 at Oh UT (zero rotational phase) and 
then phased to the period P = 13.7744 ± 0.0002 hr (SJ04). The maximum light-travel time 
difference between any two measurements included in the 2003 and 2010 lightcurves is less 
than 1.5 minutes (0.002 rotations) so no delay correction was required. More importantly, 
the relative uncertainty in the rotational period of 2001 QG298 translates to an uncertainty of 
0.06 rotations between the 2003 and 2010 data. The latter implies that the phase alignment 
in Figure |2] is accurate to ±0.06 in rotational phase. 

The vertical axis in Figure [2] indicates magnitudes relative to the maximum apparent 
brightness of 2001 QG298- The vertical alignment of the 2003 and 2010 data is artificially set 
to facilitate comparison of the photometric ranges in 2003 and 2010 (Am2003 and A??i2oio)- In 



our data, we find that at maximum brightness 2001 QG298 has magnitude rrir = 21.63 ± 0.02 
mag. This measurement is absolutely calibrated using SDSS catalog stars as described in 
Section 2 and the error is dominated by uncertainties in the r magnitudes of the calibration 
stars (Table [2]). SJ04 found that at maximum brightness 2001 QG298 has an absolute mag- 
nitude rriR^l, 1, 0) = 6.28 ± 0.02 mag. The expected apparent magnitude at the time of our 
observations is thus 

TJiRiR, A, a) = mnil, 1,0) + 5logRA + f3a = 21.46 ± 0.04 mag, 

where R, A and a are taken from Table [11 and /3 = 0.15 ±0.01 mag/° is the phase coefficient 
for KBOs ( Sheppard fc Jewi"tt||2002 ). To compare the expected m^ magnitude with our new 
nir measurement we used two approaches. Firstly, we converted our SDSS r magnitude to 



a Cousins R magnitude using the transformation equations in Smith et al. (2002) and the 
broadband colors of 2001 QG298 {B - V = 1.00 ± 0.04 mag and V - R = 0.60 ± 0.02 mag; 
see SJ04): we obtain an equivalent R magnitude tjir = 21.36 ± 0.04 mag. Secondly, we took 
the transmission profiles of the two filters (available from the Isaac Newton Group online 
filter database) and convolved each with the spectrum of 2001 QG298 to find an integrated 
fiux ratio 1.215. This corresponds to magnitude difference rrir — mR = 0.21 mag between 
the filters. Subtracting that difference from our rrir measurement we obtain an equivalent R 
magnitude m^i = 21.42 ± 0.04 mag. We assumed the spectrum of 2001 QG298 to be linear 
with a slope 22.6% (lOOOA)"^, as calculated from the object's B — V color. Although both 
approaches are uncertain, they agree at the l-cr level and fall within 0.1 mag of the expected 
apparent magnitude based on the object's absolute magnitude and geometric circumstances. 



3.2. The obliquity of 2001 QG 



298 



Between the time of the SJ04 observations in 2003 and the time of our 2010 observa- 
tions 2001 QG298 has travelled 16° along its heliocentric orbit. A diagram of the geometric 
circumstances of the two sets of observations is shown in Figure [3j The change in observing 
geometry between 2003 and 2010 has predictable effects on the lightcurve of 2001 QG298- 
Given the extreme variability displayed in 2003, 2001 QG298 was probably observed equator- 
on, i.e. at aspect angle (measured between the line of sight and the spin pole) ^2003 ~ 90° 
(Figure [3]). This assertion is backed by detailed models of the 2003 hghtcurve based on 



figures of hydrostatic equilibrium (Takahashi & Ip 2004 Lacerda & Jewitt 2007 Gnat & 



Sari 2010). Those simulations consider close and contact binaries in which the components 



are homogeneous and strengthless, and orbit each other in tidally locked, circular orbits. 
Under these conditions, the binary components assume shapes that balance gravitational 



(including mutual tidal) and rotational accelerations. Furthermore, the exact shapes of the 
binary components, their mutual orbit period (= spin period), and their density are all inter- 
dependent. The practical result is that the density can be inferred from the spin period and 
the component shapes. The simulations search the space of hydrostatic equilibrium binary 
solutions and orientations (aspect angles) for the combination that best fits the observed 
lightcurve and then use the spin period to extract the density. This technique has been 
tested on the observationally well-characterized Jovian Trojan contact binary (624) Hektor 



(Lacerda & Jewitt 2007), using observations at multiple geometries along its 12-yr orbit 



(Dunlapfc Gehrels 1969) 



In the case of 2001 QG298, all attempts to simulate the SJ04 lightcurve yield equivalent 
results and agree upon a nearly equator-on geometry, and a low bulk density p ~ 0.6 — 0.7 
g/cm~^ ( Takahashi fc Ip||2004 : Lacerda fc Jewitt||2007 Gnat fc Sari||2010 ). For our purposes. 



we will use the models of 2001 QG298 found by Lacerda & Jewitt (2007 hereafter LJ07); 
the SJ04 data and respective LJ07 model lightcurves are shown in Figure [Tj The binary 
models are rendered in Figure |4} LJ07 make use of the Roche approximation which calcu- 
lates the deformations on each component independently, assuming the other component's 



mass is concentrated on a point ( Chandrasekhar 1963 Leone et al. 1984). In the Roche 



binary approximation the binary is composed of two triaxial ellipsoids described by their 
axis ratios, their mass ratio, and their orbital separation. LJ07 include the effects of limb 
darkening when simulating the Roche binary lightcurve. They consider two extreme light 
scattering laws: a lunar-type (Lommel-Seeliger or backscatter reflection) law, appropriate for 
low albedo, rocky surfaces, and an icy-type (Lambertian or diffuse reflection) law, suitable 
for high albedo surfaces. The lunar law leads to negligible limb darkening and so generates 
symmetric lightcurves for which both mutual eclipses produce the same signature. The icy 
law causes considerable limb darkening and so produces asymmetric lightcurves if the binary 
components have different masses/sizes. The SJ04 hints at a slight difference between the 
two lightcurve minima (Figure [I]), intermediate to the extreme lunar and icy cases, so LJ07 
identify two models, one for each scattering law, that bracket the data (see Figure [I]). The 
model parameters are listed in Table |3j 

We want to investigate how the lightcurve is predicted to change between 2003 and 2010 
depending on the binary's obliquity, e (see Figure [3]). The obliquity is the angle between the 
outer (heliocentric) orbit and inner (mutual) orbit planes. We shall refer to the normals to 
the outer orbit and inner orbit planes of the binary as orbit pole and spin pole, respectively. 
Assuming that 2001 QG298 was observed exactly equator-on in 2003, then the spin pole must 
lie on the plane perpendicular to the 2003 line of sight (aspect angle 6^2003 = 90°) and the 
obliquity is simply the sky-projected angle between spin pole and the normal to the ecliptic. 
[Note: The Earth is less than (1 AU/31 AU) (ISOVvr) ^ 1.8° away from the orbit plane of 



2001 QG298 as seen from the KBO, which is neghgible.] In this configuration, the hghtcurve 
has the same shape independent of the obhquity. In 2010 the spin pole of 2001 QG298 makes 
an angle with the line of sight ^2010 = arccos(— sine sinAi/) where At/ = 16° is the angle 
through which 2001 QG298 has moved along its orbit since 2003. The 2010 aspect angle can 
take values in the range 90° < ^2010 < 90° + 16° depending on its obliquitjQ In other words, 
the equator of 2001 QG298 lies at most 16° from the line of sight in 2010. 

We took the two models in Table [3] and rendered them as they would appear from Earth 
in 2010 (Figure |4]) following the procedure in LJ07. The resulting 2010 model lightcurves 
are shown in Figure [5j As expected, at obliquity e = 0° the lightcurve displays no change 
in photometric range, Am,2oo3 = ^^^2010- However, the lightcurve phase is shifted (delayed) 
by 16°/360° ~ 0.044 as the object must rotate an extra 16° in order to appear at the same 
rotational phase as seen from Earth. As the obliquity increases, the 2010 photometric range 
decreases, down to a minimum of A?7i2oio ~ 0.7 mag for e = 90°. The rotational phase 
shift also decreases as e increases and for e = 90° the shift is zero and the 2003 and 2010 
lightcurves appear in phase. As the obliquity approaches e = 90° the 2010 lightcurve looks 
nearly identical for the lunar- and icy-type surfaces. This occurs because the minimum 
cross-section configuration becomes less affected by limb darkening effects (see Figure |4]). 

Figure |6] shows how the observed 2003 and 2010 photometric ranges for 2001 QG298 com- 
pare with the model predictions. Figure [7] overplots our 2010 data on the model lightcurves 
from Figure [5j The first point to notice is that the new data are consistent with and support 
the idea advanced by SJ04 that 2001 QG298 is a contact binary (or at least a very elongated 
object), now being observed at a slightly more pole-on geometry than in 2003. Secondly, 
we find that the new data are inconsistent with a low obliquity. In fact, the data points 
with smaller error bars are best explained by an obliquity e = 90°, both in rotational phase 
and photometric range (Figured). We find the lowest reduced x^ value (x^ = 1.7) for the 
e = 90° model. Since a unit increase from the minimum reduced x^ roughly brackets a 



1-cr confidence interval (Press et al. 1992), our data suggest an uncertainty in the obliquity 
smaller than Ae = ±30° (see Table 111). This estimate neglects the effect of the uncertainty 
in the spin period and rotational phase mentioned in Section 3.1. To address this issue, 
we used a Monte Carlo simulation that combines the photometric and timing uncertainties. 
The goal is to quantify the probability that our data are best explained by a given obliquity 
range. We generated 1000 replicas of the 2010 lightcurve and shifted each by a random phase 



^Note that we are considering the case (shown as thick soHd vectors in Figure k3| in which the pole vector 
moves away from the observer between 2003 and 2010; the opposite case in which the pole moves towards 
(and 90° — 16° < ^2010 ^ 90°) the observer (thin dashed vectors in Figure Isl) produces a similar effect on 
the 2010 lightcurve. 



drawn from a normal distribution with zero mean and standard deviation 0.06 (see Section 
3.1). Then, each rephcated data point was nudged independently in the vertical (magni- 
tude) direction by a random amount taken from a normal distribution with zero mean and 
standard deviation equal to its error bar. We found that 58% of the artificial lightcurves 
are best explained (in terms of x^) by the model with e = 90°. These are followed by the 
e = 75° model with 12%. We conclude that 70% of the artificial lightcurves sampling the 
uncertainty region are best explained by models with obliquities in the range e = 90° ± 30°. 
Finally, Figure |8] shows that for e = 90° both surface types considered in LJ07 produce simi- 
lar results implying that the 2010 data are not adequate for discriminating between the two. 
We conclude that the 2003 and 2010 lightcurves are simultaneously explained by a single 
contact binary model with large obliquity, e = 90° ± 30°. The simultaneous fit is shown in 



Figure 10 



Discussion 



The large obliquity of contact binary KBO 2001 QG298 is not unusual; it is a property 
that is shared by other extreme shape objects in the solar system. Table [5] lists the subset 
of bodies from the Planetary Data System asteroid lightcurve database ( Warner et al.|[2009 ) 
with maximum photometric variability larger than Am = 0.9 mag (indicative of extreme 
shapes) and effective diameter larger than De = 100 km. We have ignored smaller objects 
as they were probably more strongly affected by collisions and non-gravitational effects that 
could interfere with their spins. The bodies in Table [5] are well studied and known to 
have extreme shapes, and they all have obliquities in the range 50° < e < 180° — 50°. 
Although suggestive, the obliquities of the objects in Table [5] are nonetheless consistent with 
an isotropic distribution. Randomly oriented spin poles fall in that same range of obliquities 
with probability p = cos 50° = 0.64. Hence, the chance that all four objects in Table |5] 
fall in that range, p^ = 0.17, remains high. The sample size is too small to draw definite 
conclusions. 

The origin of contact binaries remains unexplored. KBO binary formation theories have 



focussed on the more easily identified resolved pairs (Weidenschilling 2002 Goldreich et al. 



2002] IFunato et ar]|2004t [Astakhov et~aL][2005| [Schhchting fc Sari|[2008al [Nesvorny et~aL 



2010 ). The general idea behind wide-binary formation models is to identify a mechanism that 



removes energy from a close but unbound pair of KBOs, converting it into a gravitationally 
bound binary. Contact binaries are perhaps an end state of initially "hard" (binding energy 
of the binary larger than the mean kinetic energy of single objects) wide-binaries that have 



been progressively further hardened by close encounters with passing, single bodies (Heggie 



1975 Goldreich et al. 2002). Recently, Perets & Naoz (2009) proposed a contact binary 



formation mechanism based on the combined effects of Kozai oscillations and tidal friction. 



Kozai oscillations (Kozai 1962 Fabrycky & Tremaine 2007) occur in hierarchical triple 



systems (the Sun and the KBO binary in this case) with large mutual orbital inclinations. 
In the solar system, the mutual inclination between the inner and outer orbits of a binary is 
by definition the obliquity, e. If the initial obliquity Eq is larger than the Kozai critical angle, 
typically Ec ~ 40°, the binary orbital eccentricity and obliquity will oscillate in an anti- 
correlated fashion. As the obliquity falls to a minimum of e = 0° the eccentricity attains 
its maximum value which depends on Eq and approaches unity for Eq ~ 90°. Such large 
eccentricities will bring the binary components very close (possibly even leading to collisions) 
allowing tidal friction to dissipate energy and gradually shrink the binary semimajor axis. 



Perets & Naoz (2009) propose that repeated episodes of Kozai oscillations and tidal friction 



could eventually lead to very compact systems - possibly contact binaries or highly elongated 
objects - with eccentricities e ~ and inclinations 40° < Eq < 180° — 40°. 



4.1. The fraction of contact binaries revisited 



Current estimates of the intrinsic fraction of contact binaries among the Jupiter Trojans 



and the KBOs assume that their spins are isotropically oriented (LJ07; Mann et al. 2007). 
Since contact binaries are identified from their high variability only if observed at an aspect 
angle sufficiently close to 9 = 90°, the apparent fraction must be corrected for the probability 
of such favorable viewing geometry. LJ07 find minimum aspect angles for positive contact 
binary identification of 6'iunar = 81° and 6'icy = 71° for lunar- and icy-type surfaces, respec- 
tively. The probability that a randomly oriented contact binary is observed at aspect angle 
d > 6'min is given by cos6'min (Lacerda &; Luu 2003) and amounts to piunar = cos 81° = 0.15 
and Picy = cos 71° = 0.33 for the extreme surface types considered. Based on the discovery 
of 2001 QG298 in a sample of 34 KBOs (see SJ04), LJ07 calculate lower limits to the fraction 
of contact binaries of /lunar > 20% or /icy > 9% depending on the assumed scattering law. 

The estimated abundance will be higher if most contact binaries have high obliquities. 



Figure 10 illustrates the solid angle region available to spin poles with high obliquities, in 
the range 40° < e < 180° — 40°, and the fraction of that region detectable from Earth. 
We numerically integrated the "detectable" and "allowed" solid angles and found their ratio 
(which equals the detection probability) to be Piu^ar — O-^ ^^^ PLy = 0-24 for the two 
surface types. The lower detection probabilities with respect to the isotropic case imply 
larger intrinsic abundances by 0.15/0.11 = 36% (lunar scattering) and 0.33/0.24 = 38% (icy 
scattering). Therefore, if contact binaries have obliquities in the range shown in Figure 10 



10 



then the lower hmit abundance would increase to /lunar > 27%, assuming lunar-scattering, 
and /icy > 12%, assuming icy scattering. 



5. Conclusions 

We have measured the visible lightcurve of Kuiper belt object 2001 QG298 and compared 
it to similar data collected in 2003. Our findings can be summarised as follows: 

1. The 2010 lightcurve of 2001 QG298 has a peak-to-peak range of Am2oio = 0.7±0.1 mag, 
significantly lower than the photometric range in 2003, Am2oo3 = 1.14 ±0.04 mag. The 
2003 and 2010 lightcurves appear in phase when shifted by an integer number of full 
rotations. 

2. The change between the 2003 and 2010 lightcurves is most simply explained if 2001 QG298 
possesses a large obliquity, e = 90°±30°. In that case, the lightcurve photometric range 
should continue to decrease, reaching a minimum of Am ~ 0.0 — 0.1 mag in 2049. 

3. Current estimates of the fraction of contact binaries in the Kuiper belt assume that 
these objects have randomly oriented spins. If, as 2001 QG298J contact binaries tend 
to have large obliquities their abundance may be larger than previously believed. 
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Fig. 1. — KBO 2001 QG298 lightcurve data (black points) obtained in 2003 by Sheppard 



& Jewitt (2004). Fits based on Roche binaries are shown as sohd hnes (Lacerda & Jewitt 



2007). Blue (dashed) line assumes an icy-type surface and red (solid) line assumes a lunar- 
type surface. The model binaries are rendered in Fig. |4j 
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Fig. 2.— Lightcurves of KBO 2OOIQG298 in 2003 (open circles; Sheppard fc Jewitt 2004) 



and 2010 (filled symbols; this work). Red circles, green squares and blue stars indicate data 
from 2010 August 15, 16 and 17, respectively. The photometric range has decreased from 
Am2oo3 = 1-14 ± 0.04 mag in 2003 to Amsoio = 0.7 ± 0.1 mag in 2010. The 2003 and 2010 
lightcurves appear in phase. 
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Fig. 3. — Geometric circumstances and relevant angles used to describe the lightcurve ob- 
servations of 2001 QG298 in 2003 and 2010. The orbit plane of 2001 QG298 is shaded and la- 
belled. The location of the Earth (black dot) and the lines of sight (dashed lines "LOS2003" 
and "LOS2010") towards 2001 QG298 are also indicated. The positions of 2OOIQG298 in 
2003 and 2010 are marked by red boxes. A few possible spin pole orientations (obliqui- 
ties e = 0°, 30°, 60°, 90° measured from the normal to the orbit plane, as indicated) for 
2001 QG298 are illustrated as thick solid arrows emerging from the object. Thinner dashed 
lines correspond to symmetric pole orientations (w.r.t. e = 0°) which produce similar 
lightcurves to their counterparts and cannot be observationally distinguished. Under the 
assumption that the 2003 aspect angle 6^2003 = 90°, the 2010 aspect angle lies in the range 
90° < 02010 < 90° + 16° (or equivalently 90° - 16° < ^2010 < 90°) depending on the obliquity. 
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Fig. 4. — Renderings of the Roche binary models of 2001 QG298 (Lacerda & Jewitt 2007) 



assuming lunar-type (top two rows) and icy-type (bottom two rows) surface scattering. The 
models are rendered as seen from Earth in 2003 (rows 1 and 3 from top) and in 2010 (rows 
2 and 4 from top) assuming that 2OOIQG298 has obliquity e = 90°. Successive rotational 
phases are displayed from left to right to simulate rotation and to show maximum and 
minimum cross-section configurations. 
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Fig. 5. — 2001 QG298 lightcurve model predictions for 2010 (solid, colored lines) based on 
the 2003 models (dashed lines). The predicted lightcurve range and phase depend on the 
assumed obliquity e. Lunar-type (icy-type) scattering models are shown in the top (bottom) 
panel. 
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Fig. 6. — Predicted photometric range for 2001 QG298 as a function of time and obliquity. 
Measured 2003 and 2010 photometric ranges are shown as large gray points with l-cr er- 
ror bars. Table |4] lists the Am values as a function of obliquity at the time of the 2010 
measurement. 



20 



0.0- 





0.2: 


-(— > 


. 


• 1—1 




bX) 


0.4- 


c^ 


. 


S 


0.6- 


<D 


■ 


> 




• 1—1 


o.s- 


OS 






. 




1.0- 



1.2- 




0.4 0.6 

Rotational Phase 



1.0 



Fig. 7. — Comparison between the 2010 lightcurve data (black dots with error bars) and 
the model predictions shown in Fig. |5| Only the lunar-type scattering models are plotted. 
Different lines assume obliquities e = 0°, 30°, 60°, 90°. Respective reduced x^ values for the 
fits are x?cd = 18-8' 12.3,4.4, 1.7, favoring the solution e = 90°. 
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Fig. 8. — Negligible difference between the lunar- (red, solid) and icy-type (blue, dashed) 
scattering model lightcurves. 
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Fig. 9. — Lightcurves from 2003 and 2010 (as shown in Fig. pi) simultaneously fit by the 
model in Fig. |4] (top two rows). The model assumes obliquity e = 90°. 
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Allowed 




Fig. 10. — Solid angles used to calculate the probability that a contact binary with obliquity 
in the range Ec < s < 180° — Ec {Ec = 40°) will be detected from its high photometric 
variability. The line of sight is indicated by a thick, black arrow labelled "LOS". Solid angle 
regions of interest are marked on the surface of a sphere surrounding the contact binary 
(marked by a black dot). The solid angle region where spin poles are allowed with uniform 
probability is shaded in light gray. That region is defined as the full sphere minus polar caps 
extending an angle Ec from the either pole. Darker gray regions mark spin pole orientations 
that are detectable from Earth. Those regions lie more than ^min from the LOS. The ratio 
of these detectable and allowed solid angles equals the detection probability. 
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Table 1. Journal of Observations. 



UT Date i?" A^ a" Filt. '^ Seeing *= Exp. -^ Notes 

[AU] [AU] n ["] ¥\ 



r 


0.9-1.5 


600 


photometric 


r 


1.3-1.7 


600 


clouds 


r 


0.9-1.3 


600 


thin cirrus 



2010 Aug 15 31.7618 31.1064 1.41 
2010 Aug 16 31.7618 31.1053 1.41 
2010 Aug 17 31.7618 31.1042 1.41 



Note. — All observations conducted at the 2.5-ni Isaac Newton Telescope. 

'^Heliocentric distance in AU; 

^Geocentric distance in AU; 

'^Phase angle in degrees; 

'^Filters used; 

*^Typical seeing in arcseconds; 

^Typical integration time per frame in seconds. 
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Table 2. Calibration Stars. 



SDSS Name 


r" 


g_^.b 




mag' 


mag 


SDSS J005028.45+020005.7 


18.473 ±0.016 


1.39 


SDSS J005023.64+015726.9 


18.951 ±0.025 


1.22 


SDSS J005026.24+015716.2 


17.528 ± 0.024 


0.56 


SDSS J005012.65+015638.6 


19.131 ±0.026 


1.32 


SDSS J005011.82+015538.0 


17.932 ± 0.024 


0.97 



*SDSS PSF r magnitude; 
^SDSS Model g — r color; 
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Table 3. 2001 QG298 models. 



Surface " 


q' 


B/A^ 


C/A^ 


h/a'^ 


c/a'^ 


ujVi-KGpY 


d/{A + a)f 


P' 


Lunar 
Icy 


0.84 
0.44 


0.72 
0.85 


0.65 
0.77 


0.45 
0.53 


0.41 
0.49 


0.130 
0.116 


0.90 
1.09 


0.59 
0.66 



"^Surface scattering type; 

"^Mass ratio of the binary components; 

'^Axis ratios of primary; 

'^Axis ratios of the secondary; 

''Dimensionless spin frequency of triaxial binary; 

^Dimensionless binary orbital separation; 

^Bulk density of the models (in g cm~^). 
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Table 4. Model Goodness of Fit. 





Model Photometric Range ^ 




Obliquity " 


Amo3 [mag 


ArriiQ [mag 


Xred 


£ = 90° 


1.14 




0.77 


1.7 


£ = 75° 


1.14 




0.79 


2.4 


£ = 60° 


1.14 




0.83 


4.4 


£ = 45° 


1.14 




0.91 


7.7 


£ = 30° 


1.14 




1.02 


12.3 


£=15° 


1.14 




1.13 


16.9 


£ = 0° 


1.14 




1.14 


18.8 



'^Obliquity of the model; 

^Model photometric ranges in 2003 and 2010 as- 
suming lunar-type scattering; 

•^Reduced x^ value of 2010 model fit to the 2010 
lightcurve. 



Name" 
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Table 5. Objects with Extreme Shapes. 



Group ^ De^ P" Am„,a> 
[km] [hr] [mag] 



Obhquity • 



Shape ' 



Refs. 



90 Antiope 


MBA 


120 


16.5 


0.90 


53 ±2 


Close binary 


1,2 


216 Kleopatra 


MBA 


135 


5.4 


1.20 


84 ±2 


Bi-lobed 


3,4,5 


624 Hektor 


Trojan 


230 


6.9 


1.10 


98 ±5 


Contact binary 


3,4,6,7 


139775 2001 QG298 


KBO 


250 


13.8 


1.14 


90 ±30 


Contact binary 


7,8,9 



References. — (l) |Michalowski et al.] ( |200lD ; (2) |Descamps et al.| ( [2009| ; (3) |de Angelisl ( |1995D ; 
(4) [Kryszczyhska et al.| ( [2007| ; (5) [Descamps et~aL] ( [20ll| ; (6) [Dunlap fc Gehrels] ( |I969| ); (7) 
Lacerda fc Jewitt| ( |2007D ; (8) [Sheppard fc Jewitt| ( |2004D ; and (9) this work. 



'^Object designation; 

'^Object group (MBA=main-belt asteroid, Trojan= Jupiter Trojan, KBO=Kuiper belt object; 

■^Effective diameter; 

'^Approximate spin period; 

'^Maximum observed photometric range; 

^Approximate obliquity. 

^Most probable shape. 



